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A new family of exactly solvable models is introduced. These models are gener- 
alizations of the XXZ chain where the distance among spins up (cr z -basis) cannot 
be smaller or equal to t (t = 0, 1, 2, ...). The case t = recovers the standard XXZ 
chain. The coordinate Bethe ansatz is applied and the phase diagram is calcu- 
lated. Exploring the finite-size consequences of conformal invariance the critical 
exponents are evaluated exactly at the critical regions of the phase diagram. 



1 Introduction 

The anisotropic S = h Heisenberg model, or XXZ chain, is one of the most 
studied quantum spin system in statistical mechanics. Since its exact solution 
by Yang and Yanga this model has been considered to be the classic example 
of the success of the Bethe ansatz. In this paper, considering the XXZ chain 
as a prototype we introduce a new set of quantum chains which preserve the 
property of integrability. 

The XXZ Hamiltonian on an L-site chain is given by 

H(A) = + + , (1) 

i=l 

where a x ,a y , a z are spin-^ Pauli matrices, A is the anisotropy parameter and 
we assume periodic boundary conditions. If we identify, in the a^-basis, a spin 
up or spin down at a given site i as a particle or hole at this site, the XXZ 
Hamiltonian can be written iquivalently as 

L 

H(A) =-J2 ^o(c+c m + ct +1 a + An t n l+1 )P , (2) 
i+l 
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where cf and Cj are spinless creation and annihilation operators and rij = 
c^Ci is the density of particles at site i. The projector P a in (1) forbids two 
particles to occupy the same position in space (fermion like). The models 
under consideration are obtained by a generalization of the operator P and 
are given by 



H(A) = -^P t (cfc i+1 + cf +1 c i + An i ni + t+i)Pt, t = 0,l,2,... (3) 

i=l 

where Pt now forbids two particles at distance less or equal to t (0, 1,2, ...). 
The static interaction, controlled by A, acts only when two particles are at the 
closest possible positions (i, i+t+1). The case t = recovers the standard XXZ 
chain (2) or (1). In terms of Pauli matrices we can also write the equivalent 
Hamiltonian 

1 L 

H t (A) = -- PM< + i + °> V i + i + Aafal +t+1 )P t , (4) 
i=i 

where now P t projects out from the Hilbert space the states (in the <7 2 -basis) 
where two up spins (at — 1) are at distance smaller or equal to t. In coming 
from (1) to (2) or from (3) to (4) we neglect a term proportional to a magnetic 
field. This term can be trivially added in the above Hamiltonian since the 
magnetization Yli=i °f or the total number of particles n = X^=i n i are good 
quantum numbers. In terms of Pauli matrices the operator Pq = 1 and for 
t > it is non-local and can be written in the c z -basis as 



|(1 -of) + £(1 + of) J] 1(1-0^,) 
1=1 



(5) 



The Hamiltonians (2), although looking artificial, may describe interesting 
physical situations where the states with up spins are associated with larger ex- 
cluded volume than those of down spins. This is precisely the case when we con- 
sider the time fluctuations due to diffusion of molecules with size t+1 (in units 
of the lattice spacing parameter) on a lattice. Interpreting the master equation 
associated to this problem as a Schrodinger equation in Euclidean timcJ3 the 
associated Hamiltonian is given by the isotropic ferromagnetic version of (3), 
namely H t (A = +1). The anisotropic ferromagnetic version_(A > 1) is related 
to the asymmetric diffusion of the molecules of size (t + l).tl 



2 



2 The Bethe ansatz equations 

It is important to observe that similarly to the XXZ chain the ^-magnetization 
to = Ylf a ii or t ne total number of up spins n (or the number of particles 
in (2) and (3)) is a conserved quantity for the Hamiltonian (4). Consequently 
the Hilbert space can be separated into block-disjoint sectors labelled by n = 
0,1, ... ,L/(t + 1). In a given sector n the wave functions of (4) can be written 

L L+t(t+l) L 

® = / I ■■■ ^ •■■ 51 a{xi,x 2 , ■ . . ,x n )\xi,x 2 , ■ ■ ■ ,x n ) , (6) 

£1 = 1 Xi + l=0^+t+l Xn=X n — 1+t+l 

where L = L — n + 1 — t(n — 1) and a(:£i, x n ) are the amplitudes of 
the spin configurations \x\, x 2 , x n ), having spin down at the coordinates 
Xi(i = 1, 2, n). We can calculate the above amplitudes and the correspond- 
ing eigenenergy of (4) by applying the Bethe ansatz 

a{x 1 ,x 2 ,...,x n )^Y. Ak ^^y {kpiX1+ --- +kpnXn ) > (7) 
p 

where the summation is over all permutations P = (pi,P2> ■■■,Pn) and ^{k} 
as well the wave numbers {k} = k\, k 2 , k n should be calculated from the 
eigenvalue equation for (4). The solution follows the standard coordinate 
Bethe ansatzD The validity of (7) for the case n = 1 is immediate due to 
the translational invariance of (4), and in this case the energies are Ek — 
—2 cos k — y(L — 4), with k — 2np/L (p = 0, 1, 2, L — 1). The case of general 
n can be verified as in the standard XXZ chain (t = 0). The energies and 
momentum are given by 

"A ™ 
E k = -2j]cos2fc 4 --(i-4n), P = ^fc*. (8) 

i=l i=l 

The wave- numbers {ki} are obtained by solving the Bethe ansatz equations 

" 1 —0\ P ik i _l_ pKkj+h) 

ik 3 L _ , pn-iTT tt(k 7 -k l ) i + e i q-i 



i=i 



These equations at t = give us the well-known Bethe ansatz equations for 
the XXZ chain. For t > these equations in a given sector specified by n up 
spins and momentum P — ^-p (p = 0, 1, L — 1) are the same as those of 
an XXZ chain also with n up spins but lattice size LL = L — tn, momentum 

itioncl 



P' = L _. p (mod27r) and twisted boundary condit 

of +1 ±ial +1 = ^(at±ia\), a£ +1 = of , ® = ^P- (10) 
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We will now study the phase diagram of these models. Calculating the eigen- 
spectra of (4) numerically by brute-force diagonalization on small lattices, we 
see that in absence of an external magnetic field the groundstate, for finite 
chains will belong to a sector where the magnetization M = 2n — L depends 
on the lattice size L, exclusion parameter t and interaction A, as long as t > 0. 
This is distinct from the standard XXZ case (t = 0) , where the groundstate is 
always in the half-filled sector n = , with zero magnetization for all values of 
L and A < 1. In order to have a well defined thermodynamic limit for given 
values of t and A it is necessary to work with fixed values of the magnetization 
per site m=^ = ^ — 1, or equivalently for fixed values of the density of up 
spins P = j,- We can force the groundstate to belong to such a sector with 
n = pL by adding in (4) an external magnetic field —h^2f =1 erf, which will 
depend on the values of t, A and p. 

It is important to remark that even when A = we do not have a com- 
pletely free system for t > 0, since in (9) a phase connecting different values 
of kj is always present. Nevertheless the Bethe-ansatz equations can be solved 
exactly for finite lattices, and is then instructive to calculate the eigenspectra 
at A = 0. 



3 The operator content at A = 

Let us consider H t (A = 0, ft.), with a given value of the magnetic field h = h(p) 
that produces a fixed density p = j; (0 < p < j^j) of up spins. Taking the 
logarithm of both sides of (9) we obtain 

pL 

Lk j =TT(pL-l) + J2t{k j -ki) + 2-Km j , j = 1, 2 (11) 

i=i 



where rrij (j = 1,2, ...,n) are integers. The solution of these equations are 
given by the set {kj}, 



2tt tt( P L -l)-tP 

LiT^tpj™ 3 ' ~ L(l - pt) 



kj=0 + — — mj, e= y \ M '— , j = l,2,...,n, (12) 



where {mj} are distinct integers and P is the momentum associated to the 
state 

P = k 3 ■ = Y P ( mod27r ) > P^Z. (13) 

3 = 1 
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The groundstate has zero momentum and corresponds to the set {fcj}, where 

2?r . . pL - 1 pL-l ^ pL-1 
L{l-tp)~ 



*j = ttFW. j = -^,-^ + i,...,^. (14) 



Using (14) in (8) we obtain the groundstate energy of (4), in the presence of 
an external magnetic field, 

E (L, p) = 2 sin / sin (j^) ~ K*P - l ) L ■ (15) 



The magnetic field h — h(p) is to be calculated by imposing that when L 

L 



Roc — Ea ^ ,p ^ will be a minimum as a function of p, i.e. 



'' = '^ = Mt^)-tV s (t^)- (16 > 

This result also tells us that in the absence of a magnetic field the groundstate 
should belong to the sector with magnetization mo = = 2po — 1, where 
h(po) = 0. In the case t = 0, po — \ and mo = 0. For t > we should 
solve h(p ) = numerically; for example, po = 0.3008443, po — 0.2193388 
and po = 0.07354512 for t = 1,2 and 10, respectively. As we can see these 
values cannot be predicted by symmetry arguments on (4) or (9). The leading 
finite-size corrections of (15) are given by 



E°(L, P ) 

L ~ e °° 6L 



TT2t + °( L ~ 2 )i ( 17 ) 



where 



- = - 2 ^ - fe) - ^ - » ■ <*-ih si " fe) ' (18) 

These relations are a clear indication that the model is critical for arbitrary 
values of t and p, since (17) is the expected behavior of a massless conformally 
invariant critical point with c and £ being the conformal anomaly and sound 
velocity, respectively.!!] The sound velocity can be calculated from the energy 
momentum dispersion relation. The lowest eigenenergy in the sector with 
n = pL spins up, having momentum ^-p (p = 1,2,...) can be obtained by 
choosing the particular solution of (11), 

2tt 2?rf 
k J = TT, 73-?' 
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L{l-tp) J (l-tp)L 
P L - 1 PL - 1 pL - 3 pL-1 

J = n ' o + !>■■■> n ' o +P' ( 19 ) 
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and the corresponding energy is given by 
E^(L,p)=E°(L,p) + 

4 S in ( *(PL+P- V ) sin ( -II—) + ( L -i) . (20) 
V (l + tp)L J \L(l-tp)J { 1 y 1 

This gives us, for large L, a linear dispersion relation with sound velocity 



l-tp \l-tp 

and consequently the conformal anomaly is c = 1 for all values of t and p. The 
conformal invariance of the infinite critical system, beyong relation (17), also 
predicts the finite-size corrections of excited states.Q For each primary operator, 
with dimension x$ and spin s^, in the operator algebra of the critical system 
there exists an infinite tower of eigenstates of the Hamiltonian, whose energy 
.Em / and momentum P^ in a periodic chain are given by 

Et„AL) = J B + ^(x, + m + m')+ O (E- 1 ), 

Pt,m> = C^+m-m')^, (22) 

where m, m! = 0, 1, 2, . . .. The groundstate corresponds to the identity opera- 
tor xi = 0, Si = and the tower of states E^~^) = 2jE (pi— m') in our model are 
obtained by decreasing (increasing) the values of the quasi-momenta fcj in the 
set {ki} given in (14). After the calculation of the eigenenergies corresponding 
to many of these descendants we convince ourselves that the degeneracy of 
such levels are given by d^~ m , = P(m)P(m') where P(m) is the number of 
integer partitions of m. Another conformal dimension in the sector n = pL 
can be obtained by taking negative (positive) values of kj in the groundstate 
root configuration (14) and replacing them by positive (negative) values. For 
example, the lowest of these eigenenergies corresponds to the solution 



2tt 

[l~p~t)L' 



k o = ~ — htfO'-p*)) 



pL-l pL-1 pL-1 
J = o — + 1 ' o — + 2, — = — + 1, (23) 



and has the value 
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and a macroscopic momentum P = 2np. Comparing the finite-size corrections 
of (24) with (22) we see that this energy corresponds to an operator with 
dimension xq.i and spin sq,i given by 

XQ.rn = (1 - pt) 2 m 2 , S ,m = nffl = pLfh , (25) 

respectively. The dimensions Xo^fn, with m positive (negative) integer are ob- 
tained by taking rh of the lowest (highest) values of kj in (14) and replacing 
them by the lower (higher) available positive (negative) values. As in the case 
of the identity operator the descendants of the above operators, with dimension 
x.m + to + to' (to, to' e Z) are obtained by increasing (decreasing) the positive 
(negative) values of {kj} in the configuration that produces Xo,m- As before 
we convince ourselves that the degeneracy of these states is P(m)P(m'). 

There exist also dimensions associated to eigenstates produced by the ad- 
dition (n = —1, —2, . . .) or subtraction (n = 1, 2, . . .) of spins up in the ground- 
state sector. In this case n = pL + n and the lowest eigenenergy is obtained 
from the set {kj} with 



kj r/1 ^ , j — Jmax? jmax + 1, . . . , jmax Ijjmaxi (^) 



2tt 

L(l ~p~t - n) 

where j max = pL+ "~ 1 . These are zero- momentum eigenstates with energies 



E n = -2sin( "( pL \"'_ I /sin 



L{l-pt)-n)' \L(l-pt)- 
— h(2(pL + n) — L) . (27) 

The finite-size corrections, when compared with (22) give us the conformal 
dimensions 

- 2 

n 

X ^ = MTW (28) 

In the general case, where we combine {kj} configurations which give (25) 
and (26) we obtain the conformal dimensions 

Xn,rh = n 2 X p + , X p = — - ^ 2 n,me Z (29) 

of primary operators with spin Sn,m = (pL + n)fh. If we calculate the eigenen- 
ergies which correspond to the descendants with dimensions x™^™ = i B lS + 
to + to' we obtain their degeneracy (i™^ 1 = P(m)P(m'). The dimensions 
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(29) are the same as the dimensions x^ fh = + nKra? of operators with 
vorticity rh and spin- wave number n in a Gaussian modelQ with coupling con- 
stant K = (1 — pt) 2 /it. The degeneracy d™'^ imply that the operator content 
of if* (A = 0) is given in terms of U(l) Kac-Moody characters. The situa- 
tion is completely analogous with that of the standard XXZ chain (t = 0) 
in its massless regime.B However, expression (29) shows different dependences 
on the density p (or magnetization), depending on the value of t. In particu- 
lar for A = the dimensions are p-independent only in the free-fermion case 
t = A = 0. The result (29) shows that the net effect of the parameter t > 
at A = is similar to that of the anisotropy parameter A ^ in the stan- 
dard XXZ chain (t = 0). In both cases the exponent X p changes continuously 
with the density p of up spins. It is interesting to remark that when p — > 
X p = 1/4, for all values of t. This independence ot t arises since the restric- 
tions imposed by t > are not important due to the small amount of up spins 
in this limit. The other limit of X p at high density, p — ► l/(t + 1), can also 
be understood as follows. If we start with the maximum density p = l/(t + 1) 
and we take out one spin up, the "hole" produced will have the same effective 
motion as that produced by adding t + 1 "particles" in the low-density limit. 
Then according to (28) we should obtain the limiting value X p = (t + l) 2 /4. 

4 The operator content at A ^ 

Let us consider the case where The Bcthe ansatz equations (9) supple- 

mented by numerical calculations give the same phase diagram for all values 
of t. For A > 1 the groundstate is it + 2)-degenerate, where (t + 1) states are 
in the sector with maximum number of up spins (n = y^j) and the last in the 
sector with no spins up (n = 0). For A < 1 we expect a massless behavior sim- 
ilar to that in the XXZ chain. In fact, exploring the similarity of our equations 
for arbitrary t, with those of the XXZ chain (t = 0), we can apply the method 
of Ref. 10 to extract the finite-size corrections for arbitrary values of A and p. 
The dimensions we obtained are given by (29), where X p is now given by 

X P = J(l-tp)-V 2 (f/o), (30) 

where Uo and rj(U) are obtained by solving some integral equations. These 
equations for — 1 < A = — cos 7 < 1 are given by 

m) . 1-S2 If" ,:MW„ | tf ,,3 1 ) 



27T cos U — cos 7 2n J_jj cosh(C7 — IP) — cos(27) 
V{U) + J_ Uo cosh(t/-^)-cos(2 7 ) ^ ' (32) 
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while for A = — cosh 7 < — 1 they are given by 

1 sinhA 1 sinh(2A)Q(C/Q 



27r cosh A — cos [/ 27r ,/_r/ cosh(2A) — cos([7 — E/ 7 ) 
/ rr x 1 f U ° smh(2AWt/') 
1 = ^ )+ 27y_ c/o cosh(2A)-cos( t /-^) d[/ - (34) 

In Eqs. (31)-(34) Uq is calculated by imposing 




q(u) = <! n*+D P (35) 

The value of X p is obtained by numerically solving Eqs. (31)-(35). In Fig. 1 
we plot X p at A = — ^ for t = 0, 1, 2 and 3, while in Fig. 2 X p is shown for 
< = 2 and some values of A (1 < A < 1). 




0.5 

P (t + 1 ) 



Figure 1: The value of X v as a function of p = ^ at A = — for t = (a), t = 1 (b), 
t = 2 (c) and t = 3 (d). 

These results coincide with those calculated earlier for the standard XXZ 
chain (t = 0)£!J At low density (p — > 0) and high density (p(t + 1) — > 1) 
the exponent X p = 1/4 and = (i + l) 2 /4, respectively. These results are 
the same as those of the A = case, since in these limits the distribution of 
"particles" at p — > or "holes" at p — > l/(t + 1) are so diluted as to render 
the effect of the interaction term A negligible. From (30)-(34) it is possible 
to relate X p (t, A) for different values of t and p. In particular at the point 
p = 1/(2 + t), which corresponds to the "half-filled" density for these models, 
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Figure 2: The value of X p as a function of p = j. for t = 2 and anisotropy — 1 < A = 
- cos 7 < 1, with (a) 7 = Jjy, (b) 7 = f , (c) 7 = f , (d) 7 = f and (e) 7 = f . 

we can show that 

*p(*. 2^*) = l S + 2 )%(* = 0,p= i) = + 2) V - 7) , (36) 

for — 1 < A = — cos 7 < 1, where we have used the known resultsi for the XXZ 
chain in the absence of an external magnetic field. 

In Fig. 3 we plot X p for t = 3 and for some values of A < — 1, while in 
Fig. 4 X p is shown at A = — 2 for t = 0, 1, 2 and 3. 

As before in the low density (p — > 0) and high density (p — > l/(i + 1)) 
limits we recover the A = case. In Fig. 3 we also notice that at the "half- 
filling" density p = 1/(2 + 1) the exponent X p = (t + 2) 2 /4, independently 
of the value of A. A physical explanation for this behavior is the following. 
Strictly speaking at the "half-filling" point p = l/(t + 2) the model has a gap 
for A < —1, as long as the magnetic field is zero. When the magnetic field 
increases, reaching a critical value h — h c (A, t) this gap vanishes. Therefore 
the main physical effect at p = 1/(2 + 1) is just the vanishing of the gap and 
not the interaction A itself. In the limit A — > —00 and p < 1/(2 + t) the 
model with a given value of t behaves effectively as a model with t' = t + 1 
and A = 0, since the spins up have infinite repulsion and will be located 
in positions with no interactions among them (distances greater than t + 1). 
This can also be verifyed directly by setting A — > —00 in the Bethe ansatz 
equations (9). Consequently for A — > —00 and p < 1/(2 + t) we should have 
X p = l/[4(l-p(i+l)) 2 ]. In particular at the "half-filling" density p = l/(2+i) 
we obtain the value X p = (f + 2) 2 /4. In Fig. 3 the limiting curve (A — > —00), 
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Figure 3: The curves X p (p) for t = 3 and (a) A = -20, (b) A = -4, (c) A = -2. The 
expected curve at A — > — oo X p = 1/[4(1 — 4p) 2 ] is also shown. 




Figure 4: The curves X p (p) at A = -2 for (a) t = 0, (b) t = 1, (c) t = 2 and (d) t = 3. 
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X p = 1/4(1 — 4(g) 2 , for t = 3, can be compared with that of A = —20, obtained 
by solving (30) and (33)-(35). For A — > — oo, but p > 1/(2 + t) this argument 
does not apply and we should recover the value X p = (t+l) 2 /4 at p = l/(i+l). 
We believe the results reported here in the A < — 1 region are unknown even 
for the standard XXZ chain (t = 0). 

5 Generalizations 

Before concludiing let us discuss an extension of our results. We only consid- 
ered the cases where t > but what should be the quantum chain that would 
correspond to t = —1? In terms of particles, like in (2) this corresponds to the 
situation where we can have an arbitrary number of particles on a site. In fact 
we were able to find an exactly integrable quantum chain in this class, namely 

L oo oo 

/MA) = -£(E E (e;^e^ + e^e;^) 

j — 1 n—0 m— n+1 

oo 

+ ^A(n-l)^; +1 }, (37) 

n=0 

where E n ' m are co-dimensional matrices located at the lattice sites with ele- 
ments E n ' m )ij = 5 n j5 m j. The Bethe ansatz equations are given by (9) with 
t = — 1. This Hamiltonian commutes with the operator 

L oo 

j = l 71=1 

which gives the total number of particles. The Hamiltonian (37) in a given 
sector with TV particles can be written in terms of spin-S* SU(2) matrices with 
S = (N — l)/2. 
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